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ABSTRACT 

Few normal galaxies have been probed using near-infrared polarimetry, even 
though it reveals magnetic fields in the cool interstellar medium better than either 
optical or radio polarimetry. Deep H band (1.6 //m) linear imaging polarimetry 
toward Taurus serendipitously included the galaxy 2MASX J04412715+2433110 
with adequate sensitivity and resolution to map polarization across nearly its full 
extent. The observations revealed the galaxy to be a steeply inclined (~ 75°) disk 
type with a diameter, encompassing 90% of the Petrosian flux, of 4.2 kpc at a 
distance of 53 Mpc. Because the sight line passes through the Taurus Molecular 
Cloud complex, the foreground polarization needed to be measured and removed. 
The foreground extinction Ay of 2.00±0.10 mag and reddening E{H — K) of 
0.125±0.009 mag were also assessed and removed, based on analysis of 2MASS, 
UKIDSS, Spitzer, and WISE photometry using the NICE, NICER, and RJCE 
methods. Corrected for the polarized foreground, the galaxy polarization values 
range from zero to 3%. The polarizations are dominated by a disk-parallel mag- 
netic field geometry, especially to the northeast, while either a vertical field or 
single scattering of bulge light produces disk-normal polarizations to the south- 
west. The multi-kpc coherence of the magnetic field revealed by the infrared 
polarimetry is in close agreement with short wavelength radio synchrotron ob- 
servations of edge-on galaxies, indicating that both cool and warm interstellar 
media of disk galaxies may be threaded by common magnetic fields. 

Subject headings: Galaxies: magnetic fields - Infrared: galaxies - ISM: magnetic fields 
- magnetic fields - polarization - techniques: polarimetry 
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Introduction 



What roles do magnetic fields play in the interstellar media (ISM) of galaxies, and 
how are those fields created and sustained? How common are galaxy-wide magnetic fields 
in normal galaxies like the Milky Way? Radio observations have revealed a wealth of 
information about magnetic fields by tracing total and polarized syn chrotron ernission and 



Farad ay rotation in the warm, diffuse gas in galaxies (see reviews by 



Beck et al. 



1996 



Beck 



20091 ). but there are few observational techniques 



cooler, star-forming ISM in normal galaxies (e.g. 



capab 



Jones 



e of probing magnetic field s in the 



1997 



Packham et al 



201l|) 



The radio observations generally show disk-pa r allel (toroida 



edge-on disk galaxies ( e.g. 



2006 



Soida et al. 



Sukumar fc Allen 



1991 



Dumke et al. 



magnetic fields for 



1995 



Krause et al 



20111 ). with transitions to more vertical (poloidal) fields in the halos. 



Optical wa velength polarizations , originally thought to arise from dichroic (selective) 



extinction ( Serkowski et al. 



Scarrott fc Draper 



Fendt et al 



(119961 ) 



19751 ) that traced magnetic fields (e.g., 
19961 ). show polarization orientations perpendicular to galaxy disks, 
carefully compared their optical polarimetry to radio synchrotron polarimetry for three 
highly-inclined disk galaxies and concluded that single scattering dominates at optical 



wavelengths. As a result, magnetic information is not revealed by optical 



of such galaxies. In the near- infrared (NIR), scattering is less important (jWood 



polarimetry 



1997 



Wood fc Jones 



19971 ). enabling NIR polarimetry to reveal magnetic fields in the cool. 



star-formin g ISM that radio synchrotron emission and polarization often appear to avoid 



flBeck 



20091 ) 



Clemens et al 



2012al ) is 



The Galactic Plane Infrared Polarization Survey (GPIPS: 
in the process of using background NIR starlight polarization at H band (1.6 /xm) to 
probe the magnetic field of the star-forming mid-plane of the northern Milky Way. GPIPS 



reveals the field to distances in the disk of about 5-7 kpc, on arcmin scales. To put 
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GPIP S findings in context and to enable testing models of the Milky Way magnetic 



field flPavel 



2011 



Jansson fc Farrar 



2012a 



b|), NIR polarimetric observations of external, 



Milky Way analog galaxies are needed . Only a few edge-on systems have been previously 



polarimetrically detected (jjones 



19891 ) or partially mapped (jjonesl 119971 ) in the NIR. 



Recent obser yations of M51 , a fac e-on grand-design spiral previously mapped with optical 



polarir netrv dScarrott et al. 



l evels ( 



Pavel &: Clemens 



(IBingham et al 



1976 



20 



Jones 



19871 ). failed to detect H band polarization, to quite low 



]^ Infrared polari metry of active 



galaxies, such as M82 



2000 



20031 ). have revealed 



20001) and others (iRuiz et al. 
complicated, multi-origin polarization components, but the degree to which these objects 
inform our understanding of the magnetic fields in more normal galaxies like the Milky Way 
is unclear. 

A new small area, deep H band polarization s urvey of a portion of the Taurus 



Chapman et al 



2011 



) was underta ken 



20131). 



Molecular Cloud complex (distinct from the study of 

to examine relationships between magnetic fields and turbulence (iFalgarone et al. 
Serendipitously, one resolved disk galaxy was apparent in the observations. The galaxy 
proved bright enough to measure near-infrared (NIR) polarization across almost its full 
extent. Given how very few galaxy observations in the NIR meet those conditions, a more 
detailed analysis was warranted. 

However, because this galaxy (2MASX J04412715+243311, designated hereafter as 
the 'T-Galaxy' for its location in Taurus), is viewed through the gas, dust, and magnetic 



19871 ). it was 



fields associated with the Taurus Molecular Cloud complex (iHeyer et al. 
necessary to calibrate and remove the foreground Galactic effects on the polarization and 
photometric observations. This was challenging because the expected T-Galaxy H band 
polarization and the foreground Taurus polarization were likely to have similar percentage 
values. Described in the following sections are the new H band polarization observations. 
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the data reduction, the general properties determined for the T-Galaxy from the H band 
photometry, the determination of the Taurus region foreground polarization and extinction, 
the correction to the T-Galaxy polarization observations, and analysis of the nature of the 
magnetic field present in the T-Galaxy and how its field relates to those in other normal 
galaxies. 



2. Observations and Data Reduction 
2.1. Polarimetry 

Ground-based near- infrared H band (1.6 /im) imaging polarimetry obse rvations were 



conducted on the night of 2011 October 13 (UT) using the Mimir instrument (jClemens et al. 
20071 ) on the 1.83 m Perkins telescope, located outside Flagstaff, AZ. The observations 
consisted of four sets of 96 images, each of 10 s exposure time. Each set consisted of six 
sky dither positions, distributed about a hex-pattern of diameter 15 arcsec, with sixteen 
images at each dither position. Those sixteen images were obtained through sixteen distinct 
rotation angles of the cold, compound half-wave-plate (HWP) to modulate the incoming 
linear polarization stellar signals and were analyzed by a cold, fixed wire grid prior to 
detection by the 1024 x 1024 pixel InSb Aladdin III array. They provided a four-fold 
redundancy for calculating the four independent polarization position angles (IPPAs: 0, 45, 
90, 135°) used to form the total intensity normalized linear polarization Stokes U and Q 
percentage]^. The pixel size was 0.58 arcsec and the full image size was 10 x 10 arcmin. 



The details of the data coL 
described in Clemens et al. 



e ction m ethodology follow closely those employed for GPIPS, as 



fl2012ah . 



^Stokes values in this paper are formally Q/I and U/I, but are designated Q and U to 
reflect conventional usage at optical and infrared wavelengths. 
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Reduction of the Mimir data also followed the GPIPS pipeline approach (jClemens et al. 



2012al ). which included detector non-linearity correction, HWP angle-specific flat-fielding, 
and dark current corrections as well as corrections for sky transmission variations and 
instrumental polarization. The images were background-corrected in two ways. For analysis 
of the stars in the field, development of 'super-sky' corrections were used to model and 
remove the background from the images, as is normally done for GPIPS data. For analysis 
of the T-Galaxy, the same images had background zones free of emission identified and 
those pixels were used in a second order fit to model and remove the background near the 
T-Galaxy. This avoids inadvertent removal of faint, extended emission from the T-Galaxy 
that super-sky use might cause, given the small dither offsets used. All of the images were 
combined to yield a Stokes / (total intensity) image, shown in Figure [1] The T-Galaxy 
appears in the top middle of the Figure. The stars have FWHM values of 3.00 ± 0.04 pixels, 
or 1.75 arcsec. 

GPIPS stellar polarization values are normally calculated from point spread function 
(PSF) cleaned aperture photometry of the stars found in the 16 HWP images for each 
field, after registering and averaging the corresponding HWP images from all dithered 
observations. Stars in the T-Galaxy field were analyzed for their polarization properties 
using the normal GPIPS stellar processing, which yielded a catalog of polarizations, 



posit ions, and JHK photometry for the stars that matched to 2MASS (ISkrutskie et al. 



20061 ) ones. A total of 239 stars had polarization measurements completed, though less 
than 10% of those were significantly detected in polarization. This low P detection rate is 
due to the low numbers of stellar photons counted for the fainter stars and the background 
levels present in the H band. Figure [1] shows the polarization properties of fifteen field 
stars exhibiting low polarization uncertainties (apA < 15°, but also filtered such that 
cTp < 1 %). The 200-1- polarization non-detection stars are still useful, as their Stokes U 
and Q percentages are gaussian distributed, and may be averaged, with weighting, to yield 
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significant polarimetric information, as described in Section 4.1, below. 

To polarimetrically analyze the resolved T-Galaxy data, the HWP images were instead 
averaged into four IPPA images and combined to produce Stokes U and Q percentage 
and uncertainty image s. These were corrected for the Mimir instrumental polarization 



( IClemens et al. 



2012bl ). before uncertainty- weighted smoothing with a 1.75 arcsec FWHM 
gaussian (matching the seeing). The resulting U and Q images were used to calculate raw 
percentage polarization Praw and position angle P.A. images, along with uncertainties. 
The Stokes /, raw polarization percentage Praw, polarization percentage uncertainty ap, 
and polarization P.A. images are all shown in Figure |21 Each panel in the figure spans a 
35 X 35 arcsec portion of the Figure [1] field that is centered on the T-Galaxy and includes 
an unrelated, Galactic star at center left. 

The total intensity Stokes / image (Figure Oa) shows an inclined, disk-type galaxy 
with a central bulge. Along the major axis, it shows emissi on above the sky level over 



20071). The color 



about 40 arcsec, or about 10 kpc at a distance of 52.7 Mpc (iCrook et al.l 
legend in the panel indicates the range of excess surface brightness in the H band, above 
the background of 14.2 mag arcsec"^, spanning 19.6 to 15.2 mag arcsec"^ in a logarithmic 
fashion. The polarization percentage uncertainty ap image (Figure [2lc) shares the same 
shape characteristics as the Stokes / image, indicating that poisson noise dominates the 
photometric, and thereby polarimetric, uncertainties. The cxp rendering is clipped at an 
upper value of 4%, and used as the basis to generate an image mask that is applied to the 
-Praw (Fig- Elb) and P.A. (Fig. |5Jd) images. These images are displayed as colored contours 
only between the limits indicated in the color legends shown at lower right in each panel. 

The Praw panel (Figure [2lb) shows polarization detection across almost the full extent 
of the galaxy, with lower values in two zones offset from the center and higher values 
(albeit more uncertain) toward the galaxy edges. The P.A. panel (Figure [2jd) displays 
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the polarization P.A.s, which have been modified to resolve aliasing about 180°, so as to 
track behavior of the pixels across the histogram of P.A.s (not shown) more smoothly. This 
results in apparent P.A. values spanning a range of 30 to 210°, as indicated in the color 
legend in that panel. The P.A. values seem to favor polarization orientations of about 130° 
to the northeast (NE) of the photometric center and closer to 188° to the southwest (SW). 
The former is close to the elongation direction of the major axis (142°, see below). Neither 
is along the direction favored by the polarization orientations of the stars in the field 
(Figured]). Thus, although correction for the Taurus foreground polarization contribution 
will be important to generating images of P and P.A. intrinsic to the T-Galaxy, the 
apparent detections presented in Figure [2] were not produced by th.e foreground polarization. 



2.2. Spectroscopy 

NIR spectroscopy observations were obtained toward the mn = 11 mag core of the 
T-Galaxy, using Mimir and the Perkins on the night of 2012 September 25. The spectral 
range covered the full H and K bands (1.4 - 2.5/im) at a resolution of about 500, using 
a 1.2 arcsec wide by 5 arcmin tall slit. Eight spectral images with 5 min exposure time 
were dithered along the slit. Wavelength calibration used an Argon lamp. Correction 
for telluric absorption was done using observation of k Tau (a bright AOV star) through 
nearly the same airmass and sky direction. All spectra were wavelength-calibrated, using 
the 2D dispersion function obtained from the Argon images, and registered before adding 
(with outlier rejection). The k Tau spectrum was fit and modeled to remove its hydrogen 
absorption lines in order to devel op a telluric transmi ssion spectrum, using the xtellcor 



routine in the SpeXTool package (ICushing et al. 



2004j ). This was wavelength matched and 



divided into the summed T-Galaxy spectrum. The resulting atmosphere-corrected spectrum 
for the T-Galaxy is shown as Figure |3l It presents a mostly thermal, cool stellar population 



- 9 - 



appearance, peaking at about 1.65 /im. There are no strong emission or absorption lines 
seen, to a limit of 5% of the continuum. This includes a lack of Br 7 (2.166 /xm), either at 
zero velocity or at the redshift indicated for the T-Galaxy (see below). 



3. T-Galaxy General Properties 



Source Catalog (XS C: 



Jarrett et al 



This galaxy was f irst studied by ICrook et al.l (120071 ). who used the 2MASS Extended 



Huchra et al. 



2005a 



20001 ) as input to the 2MASS Redshift Survey (2MRS: 



bj) to search for groups of galaxies containing at least three members. 



The T-Galaxy was listed as the se cond rank member of the five-galaxy group designated 



number 308 by 



Crook et al 



(l2007l ). whose measured radial velocity of 3868 km s"^ places 
the T-Galaxy at a distance of 52.7 Mpc (z ~ 0.013 for if = 73 km s~^ Mpc~^). Its nearest 
group neighbor is located some 11 arcmin (170 kpc) away in sky projection, so they are 
unlikely to be interacting. 

The T-Galaxy St okes / image (F igure El a) exhibits a circular radius encompassing 90% 



of the Petrosian flux (IPetrosian 



19761 ) of 11.75 pixels, or 1.75 kpc at 53 Mpc. The similarly 



computed 50% radius lead to a Petrosian concentration index of 2.9. This 50% Petrosian 
radius is about twice the FWHM seeing, so the concentration index should be unaffected. 



Sersid f 



19631. 



19681 ) disk model 



20091 ). The galaxy 



This image of the T-Galaxy was fit with a combination 
and an exponential bulge model using the MPFIT package (IMarkwardtl I 
disk was assumed to be circular when viewed face-on to allow deprojecting the fundamental 
disk and bulge parameters simultaneously. The deprojected bulge was allowed to be 
elliptical with its own eccentricity and orientation. The line-of-sight inclination angle was 
found to be 73.7±0.1° (i.e., close to edge-on), deprojected along P.A. 141.6±0.1°, measured 
East from North. 
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The independent fits to the bulge and disk allowed calculation of the bulge-to-total 
(B/T) fiux ratio, found to be 0.24±0.04. This i s a brighter bulge t han for the Milky 



198J). The T-Galaxy disk 



Way, whose B/T ratio is in the range 0.11-0.13 ( Ivan der KruitI 
scale length was found to be 11.3±0.3 pixels (1.7 kpc) and the bulge scale length was 
0.89±0.03 pixe ls (130 pc). Milky Way di sk scale leng ths are larger, with estimates ranging 



from 2.26 kpc (Drimmel fc Spergel 



200l|) to 3.3 kpc fiLopez-Corredoira et al.ll2002h . The 



T-Galaxy bulge orientation was 142.6±5.6° (well-aligned with its disk) and the bulge 
inclination was 60.9±3.7°. The radial profile shows a small excess around a deprojected 
radius of 4 arcsec (1 kpc), which could indicate the presence of luminous spiral arms in the 
disk and confirm the spiral galaxy designation. 

The H band magnitude integrated out to the Petrosian radius is 11.35, yielding an 
absolute magnitude in this band of —22.27 prior to any extinction correction (see section 
4). The radius, concentration index, B/T ratio, and absolute magnitude are all consistent 
with the T-Galaxy being a mid-type spiral. The HK spectroscopy confirms that the bulge 
does not appear to harbor a significant active galactic nucleus (AGN). 

In all of these aspects, the T-Galaxy would seem to be a Milky Way analog, though 
smaller, of somewhat lower mass, and with a brighter bulge. 



4. Foreground Extinction and Polarization Determinations 



Because the T-Galaxy is viewed through the nearby (d = 147.6±0.6 pc: 



Lionard et al. 



20071 ) Taurus Molecular Cloud complex, which is known to exhibi 



: both significan 



extinction, in the form of dark clouds and intercloud material (e.g. 



magnetic fields, as previously tr aced through 



and NIR starlight polarimetry f lMoneti et al. 



optical (iMoneti et al. 



1984 



Goodman et al. 



Heve 



984 : 



r et al. 



1992 



1987) 



iever et al 



, and 



19871) 



Chapman et al 
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20111 ). the effects of this foreground material had to be quantified and removed to obtain 
the intrinsic T-Galaxy properties. For the foreground polarization, analysis of the Mimir 
polarimetry of the stars in the field containing the T-Galaxy was performed and evaluated. 
For the photometric extinction and contributed reddening, analysis of stellar photometry 
in a somewhat larger field centered on the T-Galaxy was performed, using several available 
data sets. 



4.1. Foreground Polarization 



Within the 10 x 10 arcmin Mimir field, 239 stars were sufficiently bright for polarimetric 
evaluation (i.e., were detected in each of the four IPPAs of the HWP). Of these, twelve 
showed good polarization signal-to-noise {P/crp > 2.5) and 227 stars showed lower 
signal-to-noise. By separating the stars into two stellar samples based on their S/N ratios, 
the effects of distance on the foregroun d polarization c ontribu tion could be examined. 



This analysis followed that described in 



Clemens et al. 



( l2012cl ). where stellar sample mean 



magnitude (which relates to S/N for poisson-limited polarimetry) was shown to be a useful 
proxy for distance. The Stokes U and Q percentages for each of these two samples were 
separately averaged, using inverse variance weighting. 

The same Stokes uncertainty variance weighting was applied to the H band magnitudes 
for each stellar sample, resulting in a weighted mean of 10.7 mag for the 12 star sample 
and 13.6 for the remainder sample. If both samples contain the same types of stars, the 
remainder sample would be on average about four times more distant than the 12 star 
sample. A high fraction of the stars in both sana ples are expected to be giants, and to be 



located kiloparsecs beyond ( IClemens et al. 



2012aJ) the much nearer Taurus Molecular Cloud 



complex. There should be relatively few stars foreground to Taurus in the Mimir field, and 
any present would likely show little to no H band polarization. 
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The Stokes Q and U percentage means for the two samples were found to agree to 
within 1.5 times their propagated uncertainties. Thus, a significant change in polarization 
across the distances probed using the stars in the field was not found. The two sets 
were therefore combined and weighted means were computed for the Stokes percentage 
parameters. The weighted mean percentage values found across all 239 stars were: 

Qfg = 0.81 ± 0.04 % ; f/pc = 0.75 ± 0.04 % ; (1) 

PpG = 1.11 ± 0.04 % ; P.ApG = 21.3 ± 1.1°, 

where the subscript 'FG' identifies the foreground contribution and Ppo has been 



Ricean-corrected to debias for the polarimetric uncertainty, following 



(1l974h . 



Wardle fc Kronberg 



These parameters best represent the foreground. Galactic polarization properties 
(localized in the dust associated with the Taurus Cloud complex) along the line of sight to 
the T-Galaxy. It should be noted that the portion of the line of sight between the farthest 
Galactic star and the T-Galaxy is not represented in our sample of stellar probes. However, 
the likelihood of that region having sufficient magnetic field tracing dust at these Galactic 
latitudes to affect the color and polarization of the T-Galaxy is low. Alternatively, the 
Galactic foreground stellar polarizations could be affected by partially canceling polarization 
components along the different pathlengths to the stars. The general coherence of the 
magnetic field orientations revealed in Figure [1] and the close agreement of the Stokes U 
and Q values for the two subsamples make this also unlikely. 

Finally, because the foreground polarization is small and the optical depth at H band 
is much less than unity, subtraction of the FG Stokes percentages from the observed U and 
Q percentage images for the T-Galaxy is sufficient to correct the polarization properties for 
the effects of the foreground Taurus layer. 

In summary, the steps taken to develop images of the intrinsic NIR polarization of the 
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T-Galaxy begin with forming the four IPPA images from the 16 HWP images and using 
these to generate Stokes U and Q percentage images and their uncertainty images. These 
were each smoothed, using a gaussian of 1.75 arcsec FWHM with inverse variance squared 
weighting from the uncertainty images. From the smoothed U and Q images, the foreground 
contributions Ufg and Qfg were subtracted. The new U and Q images and their uncertainty 
images were combined to form Pfg-raw and P. A. cor images (and uncertainties), where 
the subscript 'COR' denotes correction for the foreground contributions. Finally, the 
polarization percentage uncertainty image was used to Ricean-correct the polarization 
percentage image to yield a Pqor image. To allow comparison with radio synchrotron maps, 
the Stokes / image was multiplied by the Poor image to produce a polarized intensity (PI) 
image. All of these quantities are displayed for the central 22 x 22 arcsec of the T-Galaxy 
in Figure m and described in detail in section 5.2, below. 



4.2. Foreground Extinction and Reddening 

In addition to correcting for the foreground polarization effects, correcting for the 
reddening by the Taurus layer is necessary in order to obtain accurate colors for the 
T-Galaxy. In particular, reddening correction in the {H — K) color is needed to allow 
comparing the T-Galaxy polarim etry and du st column densities to similar values for the 



Milky Way and edge-on galaxies (IJones 



19971 ). Such an assessment can reveal whether the 



polarization measured across the T-Galaxy is dominated by dichroic (selective) extinction, 
and so able to trace magnetic fields, or whether scattering dominates and is unable to reveal 
magnetic field properties. 

Several photometric data sets were available for estimating the extinction and 
reddening along this sight line, up through the DR7 data release from t he Galactic 
Cluster Survey (GCS) of the UKIRT Infrared Deep Sky Survey (UKIDSS: [Lawrence et al. 
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20071 ). Methods for estimating reddening a nd extinction from broadband photometric data 



include Near-Infrared C olor Excess (NICE 



Lada et al 



Lombardi fc Alves 



mi), NIC 



200l[ ). and Rayleigh- Jeans Color Excess (RJCE: 



:-Revisited (NIC ER 



Majewski et al. 



20111 ). 



The Appendix summarizes the steps, tests, and potential pitfalls associated with 
estimating the foreground extinction along this line of sight from these multiple data sets 
and methods. The values of Ay and E{H — K) finally adopted were 2.00 ± 0.10 mag and 
0.125 ± 0.009 mag, respectively. 



5. Analysis 

5.1. 2MASS Color and Excess 

The Figure [T] field was observed here only through a single filter {H) using Mimir. 
So, to examine the color variation across the T-Galaxy, 2MASS images in J, if, and K 
bands were examined, and an {H — K) color map was formed. This observed color map 
was corrected for the foreground color (0.125 mag) and a second correctio n of 0.15 mag wa s 



subtracted. This latter value is that typically used in the NICE method fILada et al. 



199J) 



as representing the average intrinsic color for field stars (in the absence of a reference field 
color distribution). The net effect was the production of a color excess E{H — K) map to 
use as a proxy for dust column density across the T-Galaxy. 

The distribution of this 2MASS-traced color excess is shown in the upper left panel 
of Figure with its uncertainty shown in the upper right panel. The E{H — K) panel 
reveals weak, but significant color shifts across the face of the T-Galaxy. Two symmetric 
arm-like features are seen to the East and West of the galaxy center as regions of increased 
reddening, surrounded by regions of decreased reddening. Using the orientation and 
inclination properties described in Section 3, a deprojected and rotated version of the 
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E(H — K) map was calculated and is presented in the lower left panel of Figure [51 Axes 
for this panel are presented as linear offsets along the directions of the (projected) major 
and minor axes, based on the distance to the T-galaxy. If the two reddened arm-like 
features are spiral arms, they are located only 1-2 kpc from the galaxy center. They likely 
correspond to the photometric feature noted at four arcsec radial offset in Section 3 as 
possibly indicating the presence of spiral arms. Note that this representation only spans 
the inner 5 kpc where color uncertainties are small; the full photometric extent is about 
twice as large. 

Finally, the bottom right panel of Figure O displays a classification of zones in the 
T- Galaxy, selected for their color excess uncertainty, polarization percentage uncertainty, 
and, for two of the zones, their color excess values. The selection criteria are listed in 
the second through fourth columns of Table [U Classification names of 'Core,' 'Bulge,' 
'Interarm,' 'Arms,' and 'Outer' reflect their locations and color properties. The first two 
zones have the smallest uncertainties of color excess and polarization, but no color value 
selection. These criteria map to the brightest pixels in the T-Galaxy Stokes / image. The 
next two zones have identical color excess and polarization uncertainty criteria, but were 
split into two groups about E{H — K) = 0.17 mag. This produces an 'Arms' zone that 
encompasses at least the inner, more certain color excess and polarization locations showing 
the redder extinctions identified in Figure |5la, as distinct from the 'Interarm' zones that 
exhibit bluer colors. Finally, the 'Outer' zone collects the remaining locations for which 
the polarization uncertainty is under 1.2%. Analysis of the T-Galaxy classified into these 
five zones with the foreground-corrected polarization properties (see next section), enabled 
testing for correlations between polarization and color excess that may relate to zone 
classification and galaxy structure. 
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5.2. T-Galaxy Polarization 

The foreground corrected polarization properties for the T-Galaxy are shown for the 
central 22 x 22 arcsec in Figure HI The upper left panel (Figure |Ha) shows contours of 
polarized intensity PI overlaid on a false color logarithmic representation of the Stokes / 
image. There, the brightest peak in PI is offset to the North of the photometric center, 
while a weaker secondary peak is offset to the South. An alternate explanation is that 
the photometric peak is bisected by a depolarization feature that is distributed somewhat 
linearly and mostly parallel to the major axis, as though the bright bulge were shadowed, 
at least polarimetrically, by the intervening disk. Superposition of disk-parallel and 
disk-normal (vertical) magnetic fields could also account for the apparent depolarization. 
The PI image shows extensions along the major axis to the NW and SE, with hints of 
localized maxima at the ends of the major axis. 

The polarization percentage image (Figure IHb) reveals polarizations of about 1-2% 
across most of the extent of the T-Galaxy. There is a prominent zone of nearly zero 
polarization percentage surrounding the isolated 1% polarized feature centered at an offset 
of (3, —3) arcsec. This same feature shows distinctly different P. A. values (Figure Hlc) than 
the remainder of the T-Galaxy. Most of the T-Galaxy shows a relatively uniform P.A. of 
about 130°, which is similar to the major axis position angle on the sky. 

When the corrected polarization percentage and position angle data are combined with 
the photometric image, the result, shown in Figure Hid, reveals coherent behavior across 
many kpc of separation. The vectors shown represent the Pcor and P.A. cor behavior 
averaged over 3x3 pixels (~FWHM seeing) and rendered only every three pixels in 
each direction to preserve independence. A predominance of disk-parallel polarization is 
seen, though in the region SW of the center the polarization vectors become disk-normal. 
This SW feature may be evidence for a vertical (poloidal) magnetic field or it could 
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signify scattering of bulge light by disk or halo dust. The absence of a symmetric zone of 
disk-normal polarizations to the NE of the center would seem to favor the latter explanation. 

Overall, the NIR polarizations look remarkably similar to many of the radio synchrotron 



polari z ation maps of t 



(1l99lh 



Dumke et al. 



le dis k portions of edge-on galax ies, such as those i n 



(119951), 



Krause et al. 



fcooeh . and 



Soida et al. 



Sukumar &: Allen 



(1201 ll ). In particular, th e 



combi ned Effelsberg 100 m and VLA radio polarization map of NGC5775 by 



Soida et al. 



( 120111 ) exhibits virtually all of the same magnetic field characteristics as seen in the NIR 
for the T-Galaxy. Both galaxies show nearly edge-on disks domiriated b y disk-parallel 



magnetic fields. The NGC5775 radio map (Figure 7 of 



Soida et al. 



( 120111 )) shows the same 



disk-normal component offset from the center to the SW as seen here in Figure Hid. 

The NGC5775 radio observations are better able to probe fields in the warm ionized 
plasma and rev eal an 'X-shaped' halo field off of the disk plane (a model is sketched 



in Figure 12 of 



Soida et al 



(|201l[ )). The NIR of Figure jUd may be revealing a similar 
X-shaped halo field to the SW, but some caution must be applied, as the NIR polarization 
is more strongly affected by scattering than would apply to the radio. It is fortuitous that 
the VLA observations of NGC5775 and the NIR data reported here for the T-Galaxy have 
about the same effective resolutions (about 10-15 independent samples along the luminous 
disks), though NGC5775 is at half the distance of the T-Galaxy and is twice as large, 
physically. NIR polarimetry of NGC5775 could help establish whether the radio X-shaped 
halo fields are capable of being revealed in the NIR. 

The NIR better traces the cool, star-forming ISM in the T-Galaxy, revealing a coherent 
magnetic field with a disk-parallel orientation over extents as large as 4 kpc. In addition, 
there is weak indication of a rise in polarization percentage and P.A. flaring near the ends 
of the major axis that may hint at connections to a more vertical orientation for a halo 
magnetic field of the T-Galaxy or a decrease in the disk-parallel field strength. At present. 
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there are no radio observations of the T-Galaxy which would permit direct comparison 
with the NIR. Based on the NIR observations, a disk-parallel orientation for the B-field is 
predicted for synchrotron observations of the T-Galaxy disk at short cm wavelengths. 



5.3. Polarization versus Extinction 



In the disk of the Milky Way, sight lines to reddened stars show infrared polarization, 
with the correlation between reddening and polarization well-modeled using a mixture of 
unifo rm and random m agnetic field components distributed into multiple decorrelation 



cells ( 1 Jones et al. 



1992 



hereafter JKD). The polarization arises from the interaction of 
the magnetic field with spinning, aligned dust grains, yielding dichroic extinction of the 
background starlight. 



In the galaxy NGC 4565. lJonesl (119971 . in his Figure 4) found that points sampled mostly 
across the edge-on dust lane exhibited similar polarization and color excess behavior to the 
JKD model, while the similarly edge-on galaxy NGC 891 showed weaker polarizations. 

The T-Galaxy is less edge-on than either NGC 4565 or NGC 891. Thus the sight 
lines through the T-Galaxy disk will be geometrically shorter, though they could still 
achieve significant optical depths. This is true even at H band, because although the dust 



extinction per unit gas c olumn is much 



more than twice as high (IWood &: Jones 



ess th an in the optical, scattering efficiency can be 



19971 ). 



The five zone classifications of Figure Od were applied to the T-Galaxy (foreground- 
corrected) images of polarization percentages and color excesses to identify pixel values 
in eac h zone . There were averaged and plotted, along with summaries abstracted from 
JonesI Jl997n for NGC 4565 and NGC 891, in Figured Each of the five T-Galaxy zone 
characterizations are shown as colored crossed error bars representing the 1 a dispersions 
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of pixel values of polarization percentage and color excess for each zone. These zone 
average values and their uncertainties a re also l isted in the last two columns of Table [TJ An 
abstraction of the JKD curve, based on I Joned (119971 ). is also shown, surrounded by a gray 



band depicting the approximate 1 a range of stellar values seen in Figure 6 of JKD. 

In addition to the T-Galaxy zone averages and dispersions, values of the polarizations 
and color excesses for stars appearing in the Mimir field of view are shown in Figure |6l 
Seven of the high signal-to-noise polarization stars are represented as black open diamonds. 
The remaining, low significance polarization stars were grouped into similar NICER-traced 
extinction bins and averaged (using Stokes U and Q percentages, with inverse variance 
weighting). These are plotted as black open triangles with vertical error bars. All of the 
stellar values reflect the reddening and polarization induced on the distant stars by the 
foreground Taurus Molecular Cloud complex. 

The T-Galaxy zone averages and sample extents all fall either within or very close to 
the region populated by NGC 4565 measurements and also mostly fall on or very near the 
JKD curve. This high degree of agreement of the T-Galaxy zones with the JKD model 
strengthens the argument for dichroic extinction and a magnetic field origin for the NIR 
polarizations. 

The mean zone polarizations and color excesses are plotted against zone designation 
in Figure [71 The mean polarizations show a weak rise from the Core to the Outer zone 
of the T-Galaxy. The reddest regions are the Arms, Core, and Bulge. Though the color 
excess values could be due mostly to the zone definitions, the rise in polarization is not. 
A similar weak rise in polarization with distance from a bright central region was noted 



by iJoned (jl997l ) for the case of the nuclear region of IC 694 in Arp 299. He argued that 
this was a result of the change from poloidal to toroidal magnetic field dominance in the 
polarizing path length. Alternatively, the central region of the T-Galaxy could be suffering 



-20- 



either depolarizing due to long pathlengths through the intervening disk or a transition 
from dichroism in the NE disk to single scattering in the SW portion. 



5.4. Polarization Directions Within the T-Galaxy 



In the modehng of IWood fc JonesI (119971 ). such 'polarization nulls' near, but offset 
from photometric centers arise from an admixture of scattering and magnetic field induced 
dichroic extinction producing a 'crossed-polaroid' effect. This is most prominent for light 
originating from the far side of steeply inclined galaxies. If these effects hold for the 
T-Galaxy, then the broken symmetry revealed in the NIR polarimetric observations may 
lift the degeneracy on inclination angle and identify the near and far sides of the T-Galaxy 
disk. 



The T-Galaxy foreground-corrected polarization P.A. 
quant itatively to permit assessing the applicability of the 



map o 



Figu re Hl c was examinee 



Wood! f ll997h and 



Wood fc Jones 



( 119971 ) models. Pixels in the P.A. map were apportioned into 25 equal angle azimuthal 
wedges (relative to the photometric center) and were also subdivided into the five zone 
classifications. 

Figure |8] displays how the resulting polarization P.A. and location are related in the 
T-Galaxy. Because the Figure Hlc map was not deprojected in this analysis, the number 
of pixels apportioned into each azimuthal wedge varied. This is shown in Figure [S] as the 
lower, solid black curve and the right hand vertical axis that together report the number 
of pixels included in each azimuthal wedge. The horizontal axis advances North (0°), 
counterclockwise from West (270°), which is equivalent to the usual 'East from North' 
meaning with a 270° horizontal axis shift. The azimuth directions to the NW (320°) and 
SE (140°) major axis extents of the T-Galaxy are seen as increases in the pixel numbers in 
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those bins. 

The colored points and connecting hnes show the mean polarization P.A. values for 
each zone in each azimuth wedge, with the exception of the Core zone. The Core spans only 
four pixels, right at the photometric center, hence it supplies little azimuthal information. 
Nevertheless, the Core shows significant polarization, at an average P.A. of 147°, shown as 
the dashed violet line and error bar. Solid lines connect plotted points for the other four 
zones, though azimuth wedges lacking particular zonal pixels were not plotted. The large 
black filled diamonds connected by thick solid black lines represent the all-region average 
for each azimuth wedge. 

There are a couple of important trends in Figure [H For the half of the T-Galaxy 
located roughly North and East of the major axis, between Azimuth Wedge P.A. values of 
270 and 110°, there is a mostly linear, negative slope swing in polarization P.A. from about 
160 to 120°. A very different swing in P.A. is present in the portion of the T-Galaxy located 
South and West of the major axis, where the polarization P.A. values jump up to a mostly 
constant value near +200°. 

Two gray zones in the Figure highlight the behaviors expected for a centrosymmetric 
(scattering) pattern and for a ra dial pattern, which ca n also arise from scattering, though 



not from a single central source (IWood fc Jones 



19971 ). None of the P.A. values line up 



with the centrosymmetric zone, ruling out strong centrally-illuminated single scattering 
dominating the NIR polarization pattern for the T-Galaxy. 

The SW half of the T-Galaxy shows some of the weakest polarization percentages 
(Figure |Hb) with a region-averaged Poor of 0.4 ± 0.1%, while the NE half averages 
1.4 ± 0.1%. The zone of low SW polarization and most radial polarization P.A . orien tations 



may be evidence for the 'crossed polaroid' effect described in 



Wood fc JonesI fll997h for 



the far-side of inclined disks, though mostly affecting only a portion of the SW half. The 
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near-side of the disk, here the NE half, better permits starhght to experience magnetic field 
tracing dichroic extinction, leading to the observed polarization, while the far-side (SW) 
suffers depolarization and/or single scattering. 

There is remarkably little difference in the polarization P.A. azimuth wedge averages 
for the different zones in the T-Galaxy. The Outer zone may exhibit somewhat more 
angular scatter, but this is likely a result of the higher polarization P.A. uncertainties there. 
The strong agreement of the polarization P.A.s across the many zones and across almost 
the full extent of the T-Galaxy is strong evidence for a large-scale uniform magnetic field 
component. 



6. Discussion 



In exploiting the opportunity provided by the serendipitous NIR polarimetric 
observation of the T-Galaxy, as augmented by analyses of other existing stellar data, have 
clear understandings emerged regarding the nature of the T-Galaxy and how its properties 
might relate to those in the Milky Way and other normal galaxies? 

First, is the T-Galaxy a Milky Way analog? Is it a disk galaxy that lacks a central 
AGN (so that the polarization observed is not AGN related)? The decomposition of the 
Stokes I image into bulge and disk components argues for a significant bulge contribution 
to the total luminosity, but not one so high as to be AGN dominated. The B/T ratio 
and Petrosia n concentration index are both comparable to values seen for other disk 



galaxies (e.g., 



Benson et al. 



2007 



Lackner fc Gunn 



scale length, and bulge scale length (iMacArthur et al. 



20121). as are the surface brig htness, disk 



2003 



Ganda et al 



2009h. Optica l 



spectr oscopy of the nucleus could test for the presence of an AGN, though iCrook et al. 



(120071 ) do not comment on this. Our own HK-band spectroscopy shows a cool stellar 
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population (thermal-like) spectrum, that is unchanged when deextincted. Further, the lack 
of narrow or broad emission lines in our HK spectrum argues against the presence of a 
strong AGN. A fairly conservative upper limit of 0.15% of the T-Galaxy luminosity for 
its nuclear emission lines was arrived at through combining three times the noise in the 
spectrum with a notional linewidth of 150 km s~^. Such a low upper limit on the spectral 
line power for the nucleus of the T-Galaxy rules out normal Seyferts and other AGNs. 



The bulge to disk luminosity ratio is consistent with Sb-Sbc types f lGraham fc Worley 



20081 ). A bump at four arcsec radial offset in the photometric profile of the Mimir H band 
image used in the Sersic fits hinted at possible spiral arms. The 2MASS [H — K) color 
excess image, after foreground reddening correction, showed evidence for symmetric red 
features emanating from the vicinity of the core/bulge and extending to the T-Galaxy 
edges. These red features are likely to be tracing the dust enhancements expected for spiral 
arms in the T-Galaxy and they appear at the right radial offset to match the photometric 
bump. The 2.3 kpc circular radius to 90% of the Petrosian flux is significantly smaller than 
the extent of the Milky Way disk, but other aspects indicate that the T-Galaxy is a fair, 
though smaller. Milky Way analog. 

Second, have magnetic fields in the T-Galaxy been detected? Polarization can arise 
from many sources, though in the NIR, gen erally only sc atteri ng and dichroic extinc tion will 



dominate in normal galaxies. Modeling by 



Woodl fll997h and 



Wood fc JonesI fll997h of NIR 



polarization from steeply inclined disk galaxies explored the observational effects of each of 
scattering an d dichroic extinction alone, as well as together. When scattering dominates. 



Wood 



(119971 ) found that face-on galaxies show strong centrosymmetric polarization patterns 
and edge-on galaxies show polarizations perpendicular to the disk. Dichroic extinction tends 
to produce polarizations parallel to the disk, since the disk magnetic field they modeled was 



toroidal. The model of mixed scattering plus dichroism for an inclined galaxy, shown as 
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Figure 3 (center) in IWoodI (119971 ). contains many of the polarization features and properties 
present in Figure Hid here. 

The azimuthal variation of polarization P.A.s, shown here as Figure [H reveals both 
the negative-sloped swing of P.A. for the higher mean polarization (NE; near) side of the 
T-Galaxy disk as well as the strong positive unity-sloped P.A. swing that signifies radial 



(single scattering) for the lower mean polarization (SW; far) side of the disk. 



variation is also seen in the radio polarization map of NGC5775 ( ISoida et al 



'his azimuthal 



201l|). Similar 



quantitative characterizations in future modeling studies could help sharpen comparisons 
and interpretations. This would also address whether NIR and radio polarizations are both 
probing halo magnetic fields equally or if the NIR is affected by scattering. 

The overall match of the T-Galaxy polarization percentages and color excesses with 
the Milky Way based JKD model and with the values seen toward the edge-on galaxy NGC 
4565, shown in Figure O argue for a common magnetic field plus dichroic extinction origin 
for the polarizations and reddenings from these three galaxies. 

At present, the T-Galaxy lacks resolved kinematic observations, which are needed to 
ascertain the rotation direction. The polarization map does reveal a broken symmetry 
about the major axis, which interpreted under current models was used to predict the near 
and far-side locations of the T-Galaxy disk. The same map also showed coherent P.A. 
swings. Kinematic data would be able to test these near-side/far-side predictions. 

This galaxy was not targeted, but was contained in a field of deep Taurus stellar 
observations, and these stellar values were able to reveal and calibrate the foreground 
polarization properties. The stellar polarization properties show the orientation of the 
local Galactic magnetic field contained within the Taurus layer and could be combined 
with observations from Planck to improve the remova l of Galactic foregrounds from CMBR 



observations aimed at mapping B modes (jJones 



20031). 
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7. Summciry 

Deep H band linear polarization observations toward one 10 x 10 arcmin field in Taurus 
serendipitously contained an inclined, 53 Mpc distant, (small) Milky Way analog galaxy. 
The polarization data were sensitive enough to permit measuring the intrinsic polarization 
generated by the combination of magnetic fields and aligned dust grains across nearly the 
full extent of that galaxy. 

Because the galaxy was observed through the gas and dust in the Taurus Molecular 
Cloud complex, it was necessary to measure and remove the foreground contributions by 
Taurus to the distant galaxy polarization, extinction, and reddening. The foreground 
extinction and reddening were measured using multiple methods applied to several stellar 
photometric data sets. Corrected for this reddening, an E{H — K) color excess map of the 
galaxy revealed red spiral arms at locations in agreement with a photometric bump in the 
radial H band profile. 

Once foreground-corrected, the distant galaxy NIR polarizations were found to be 
comparable to Milky Way stellar values, indicating that they most likely are tracing the 
uniform component of the magnetic field in the galaxy disk. The field revealed shows 
remarkable similarities to the magnetic fields revealed by radio synchrotron polarimetry 
of edgc-on disk galaxies. As NIR polarimetry for this galaxy mostly arises from dichroic 
extinction in the cool, star-forming interstellar medium and radio synchrotron arises in a 
much warmer component, the correspondence of the magnetic field patterns for the two 
favors a single uniform magnetic field threading both components. 

A broken symmetry in the degree of polarization and its orientation direction with 
respect to the major axis was found, likely due to a change in the relative mix of dichroic 
absorption and scattering from the NE front and SW back sides of the disk of that galaxy. 
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The approach of using the stellar foreground polarization properties to characterize 
the local, Galactic magnetic field projected orientation could prove useful for constraining 
models aiming to correct Planck CMBR observations for Galactic contamination. 



This publication makes use of data products from: the Two Micron All Sky Survey, 
which is a joint project of the University of Massachusetts and the Infrared Processing and 
Analysis Center (IPAC) /California Institute of Technology (CalTech), funded by NASA 
and NSF; the NASA/IPAC Extragalactic Database (NED), which is operated by the Jet 
Propulsion Laboratory (JPL), CalTech, under contract with NASA; the Wide-field Infrared 
Survey Explorer, which is a joint project of the University of California, Los Angeles, 
and JPL/CalTech, funded by NASA; observations made with the Spitzer Space Telescope, 
which is oper ated by JPL/Ca l tech, under a contract with NASA; the UKIDSS project. 
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A. Foreground Extinction and Color Excess Determination 



The NASA/IPA C Extragalactic Database (NED) reported results from the model of 



Schlegel et al. 



(119981), based on COBE /DIRBE and IRAS/ISSA imag e 



using Sloan Digital Sky Survey (SDSS) data by 



Schlafiy Sz Finkbeinerl (1201 if ), for estimating 



analy sis, as updated 



foreground extin ctions toward 



about 6 arcmin (ISchlegel et al. 



:he T - Galaxy. These had an effective angular resolution of 



19981). compar a ble to the entire Mimir field of view. The 
values reported for the ISchlafiy fc Finkbeinen (120111 ) update are Ay = 4.12 7 mag and 



= 0.676 mag, the former is listed in the first row of Table [2] (the original 



1998 



Schlegel et al 



values are 



reddening, using 



isted in italics in t he se cond row of the Table). When converted to 



Indebetouw et al 



( 120051 ). the Taurus layer contribution would have been 
E{H -K) = 0.26 mag. The overall observed T-Galaxy {H - K) color from 2MASS is 0.42 
mag, as reported in NED (and was confirmed here by differential photometry in the images, 
tied to the colors of eleven 2MASS PSC stars located within two arcmin of the T-Galaxy). 
Hence, correction of the observed color by this indicated reddening would yield a galaxy 
{H — i^) color of 0.16 mag, w' hich is significantly bluer than the colors of normal spiral 



galaxies (IMannucci et al. 



200ll ). Indeed, a quick search of SDSS galaxies reveals only 3.3% 
show bluer colors (A. A. West, priv. comm.). This blue galaxy conclusion seemed unlikely, 
so several other methods and data sets were examined and the initial two rows of Table [2] 
were not used in the subsequent analyses. 

The direction to the T-Galaxy misses the most highly extincted Taurus clouds (the 
direction is near and betwee n the Eastern ends of T MC-l/HCL-2 and TMC-2/B18) showing 



the strongest CO emission (IGoldsmith et al. 



20081) in the region, and thu s is expected to be 



a low-density, low-column ISM line of sight. The iGoldsmith et al. 



( I2OO8I ) maps suggest a 



total H2 column density between about 2 — 2.5 x 10^^ cm ^ (Ay = 2 — 2.5 mag; row 3 of 
Table [2]). This extinction range is also indicated in the NICER-based map, derived from 
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2MASS JHK photometry, presented by lLombardi et al.l (I2OIOI ) and listed as the fourth row 
of Table H 

To obtain better extinction and reddening estimates, catalog data from several surveys 
were mined and analyzed for a region centered on the T-Galaxy and also for a reference 
region located 1.9° away (a[J2000] = 7 1 °125: ^[J2000] = +22 °833) The reference r egion 



was chosen, based on examination of the 



Goldsmith et al 



fl2008h and 



Lombardi et al 



fl2010h 



maps, as a direction relatively near on the sky to the T-Galaxy and showing no associated 
Taurus cloud CO emission or 2MASS extinction. Stellar photometry catalogs for both the 



T-Galaxy zone and the reference zone were searched out to 10 arcmin radii 



centers. The catalogs included 2MASS, UK 



2004( 1 archive, and the WISE flWright et al. 



DSS, the Spitzer Space Telescope (IWerner et al. 



rom those 



2OIOI ) All-Sky Data Release products. 



The resulting point source catalogs were position-matched and merged with the Mimir 
polarimetry and H band photometry several ways to test the effects of inclusion of the 
various catalogs. For exa mple, the Spitzer archive included deep InfraRed Array Camera 



(IRAC: 



Fazio et al 



20041 ) observations for a Taurus 'Legacy-Class' survey (iPadgett et al. 
20081 ) that overlapped the T-Galaxy field, but did not overlap the reference field. The WISE 
all sky survey observed both fields, but had higher uncertainties in its W2 (4.6 yum) band, 
compared to B2 (4.5 fim) for IRAC. The best photometry toward the T-Galaxy field would 
include the Spitzer data, but such a choice could introduce bias. The Spitzer IRAC/B2 
magnitudes for stars in the T-Galaxy field were compared to the matching WISE/W2 
magnitudes and found to be in excellent agreement (mean difference = 0.011 ± 0.012 mag 
for 374 stars), indicating that use of Spitzer for the T-Galaxy field and WISE for the 
reference field would not introduce meaningful bias. 

Use of UKIDSS (only the Ki band data for this region were contained in DR7) and 
Mimir deep H photometry added fainter stars beyond the 2MASS sensitivity limits, though 
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most of these did not have correspondin g Spitzer or WISE phot ometry. Given the weak 



reddening seen for Spitzer IRAC bands (llndebetouw et al 



20051 ). compared to the J, H, 



and K bands, all WISE/Wl (3.4 jum) and Spitzer/Bl (3.6 fim) magnitudes were considered 
to be equivalent L' (3.45 ^m) band values, and similarly the W2 and B2 magnitudes were 
considered to be equivalent M' (4.75 fim) band values, despite the somewhat different 
bandpasses and mean wavelengths for the bands used for the two missions. The combined 
data sets for each region were assembled using the lowest uncertainty photometric value for 
each star for each band (with sufficient departure tests from this requirement to ascertain a 
lack of important biases). 

The three main methods for obtaining extinctions and reddenings from these merged 
data sets used the data in different ways. The NICE method relies only on H and K bands, 
and infers extinction from the reddening in the {H — K) color relative either to means for 
unextincted stars, or to means of stars in an unextincted reference field. Here, the 581 
stars with {H — K) colors in the T-Galaxy region yielded an average Ay of 1.66 ± 0.12 
mag. A value of —0.14 ± 0.12 mag was found for the 475 stars in the reference field. Both 
values were obtained w hen each region's stars were compared to an average intrinsic color 



{H — K)o of 0.15 mag (ILada et al.lll994l ). The near zero extinction for the reference field 
obtained for the NICE method would appear to be confirmation of that method's utility. 
Differencing the T-Galaxy and reference zone values yielded the Ay listed in the fourth row 
of Table H 

The NICER methocj^ utilizes either all three 2MASS bands (J, H, and K), or can be 
extended to include any number of bands. Applying NICER/3 (J, H, and K bands only) 
to the reference region returned Ay = 0.71 ± 0.08 mag, in sharp contrast to the near zero 
value from NICE. Including the L' and M' band magnitudes (i.e., NICER/5) increased this 



■^Using the C. Beaumont IDL implementation 
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to 0.88 ± 0.09 mag, increasing the concern that NICE may be in error. 

The RJCE method is based on {H — M') colors, as the intr insic spread of such colors 



for unextincted stars is exceedingly small (jMajewski et al. 



20111 ). much smaller than for 



the {H — K) color used by NICE. Application of RJCE to the reference region yielded 
Ay = 0.61 ± 0.19 mag, not far from the NICER values. 

Hence, the reference region, though selected to avoid the Taurus dark and molecular 
clouds, and although showing no extinction under the NICE method, exhibits significant 
NICER and RJCE traced near-infrared color and therefore dust extinction, most probably 
associated with atomic or ionized hydrogen gas and not molecular gas. Examination of 
WISE W3 (12 /im) images revealed the presence of extended dust emission, in the form of 
IR cirrus, confirming the stellar color findings from NICER and RJCE. 



A second reference region, previously studied by iSwift fc WelchI ( 120081 ) . was selected 



to serve as the zero extinction calibrator. While still within the large Taurus region, it 
is further on the sky from the T-Galaxy than the first reference region. This may affect 
comparison of the stellar types present in this second reference region and the T-Galaxy 
field stars. Similar to the approach already described, existing photometric databases 
were mined for photometry in the five NIR bands. This new region, similar to the first 
reference region, had WISE data but no Spitzer data. The stellar colors, uncert ainties, 



and covariances needed to calibrate the NICER method (iLombardi fc Alves 



200l|) were 



obtained for the 261 stars located within 10 arcmin of the second reference center (q;[J2000] 
= 65°. 00, 5[J2000] = 17°. 00) and having low uncertainties (< 0.07 mag) for all colors. Using 
this second reference zone calibration, the T-Galaxy region stellar data set was analyzed 
using NICER/3 (i.e., JHK), NICER/5 {JHKL'W), and RJCE, with all resuhs listed in 
the final three rows of Table El 

Figure |9] presents an illustration of the application of the RJCE method. In the 
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left panel, the color-color distribution of (J — K) versus [H — M') stellar values for the 
T-Galaxy field are presented and, in the right panel, the distribution of Ay values derived 



from thes e data using R, 



show the 



Bessell fc Brett 



CE ar e shown. The solid and dashed black lines in the left panel 



(119881 ) loci for unextincted giants and dwarfs, respectively. The 
near vertical run between 0.1 to 1.2 mag of (J — K) highlights the narrow spread in intrinsic 
{H — M') colors exploited by the RJCE method. 

To compute the Ay histogram shown in Figure [9lb, the measured photometric 
uncertainties were propagated to extinction uncertainties. These uncertainties were 
typically around 1 - 1.5 mag, much larger than the bin widths typically used when 
constructing Ay histograms. To build this histogram, each Ay value and its uncertainty 
were replaced with a set of several thousand, normally distributed, Monte Carlo realizations 
and these were placed into the corresponding histogram bins. Without this probabilistic 
approach. Ay histograms become artificially noisy and subject to over-interpretation. The 
Ay histogram in Figure Ob shows that the direction to the T-Galaxy is well-represented by 
a single extinction and polarization layer, localized to the Taurus Molecular Cloud complex, 
with little to no extinction beyond Taurus to the multi-kpc distances of the Galactic stars 
probed. 

Taking into account the several methods used for analyzing the stellar data sets, 
values of Ay = 2.00 ± 0.10 mag and E{H - K) = 0.125 ± 0.009 mag (Table Wj appear to 
be the best characterizations of the line of sight to the T-Galaxy. The most conservative 
correction might consider only the stellar difference values between the T-Galaxy region 
and first reference region examined, but, as described above, that reference region (and by 
inference most of the Taurus regions mapped in CO and using NICER) seems to exhibit 
an extinction nearly as large as that found for the T-G alaxy direction. This required use 



of the second, more angularly offset reference region of 



Swift fc Welch! (|2008|). The even 
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higher extinctions found by ISchlegel et al.l (Il998[ ) and ISchlafly &: Finkbeinerl (1201 if ) through 



combining COBE, IRAS, and SDSS data (first two rows of Table [2]) may be traced to a 
strong departure of the modeled dust temperature along this direction (D. Finkbeiner, priv. 
comm.), justifying avoidance of their use here. 

For each of the NICE, NICER, and RJCE analyses, extinction maps across the 
T- Galaxy region were also developed. These did not show large gradients in extinction, nor 
particular holes or peaks at the T-Galaxy position. The maps have higher uncertainties, 
however, as only a handful of stars contribute information at each map location. Given 
the lack of new insight offered by these maps, the region mean values listed above (and in 
Table [2]) were adopted for all corrections. 
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Fig. 1. — Combined Mimir H band Stokes / image toward the field in Taurus containing 
the T-Galaxy. Linear polarizations are shown as black vectors for fifteen stars in the field 
exhibiting low polarimetric uncertainties. Vector lengths represent percentage polarization 
P, with a reference length shown in the lower left. Vector angles represent the polarization 
position angle, P. A., measured East from North in this equatorial representation in J2000 
coordinates. 
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Fig. 2. — Mimir-observed photometric and polarimetric properties of the T-Galaxy. (a. Top 
left) Stokes I image of a 35 x 35 arcsec portion of Figure [H displayed as logarithmically 
stepped false color contours of excess surface brightness [in mag arcsec"^] above sky. The 
feature at offset (16, —2) arcsec is a foreground field star (see Figure [1]). (b. Top right) 
Raw percentage polarization Praw, with contour levels stepped linearly from to 3%. (c. 
Bottom left) Uncertainty in percentage polarization crp, with contour levels of 0.25%, then 
0.5 through 4% stepped linearly, (d. Bottom right) Observed polarization position angle 
P. A., measured East from North, masked outside crp < 4%, and displayed as contours 
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Fig. 3. — Near-infrared spectrum obtained toward the bright core of the T-Galaxy. Horizon- 
tal axis shows wavelength in //m and vertical axis is in flux units, uncorrected for slit losses 
relative to the observations of k Tau. Upper traces show the spectrum obtained through the 
H and K band atmospheric windows. Lower traces indicate the uncertainties. Features at 
1.95 and 2.0//m are due to incomplete atmosphere removal. The vertical line indicates where 
the redshifted hydrogen Br-7 line should appear, if present. The lack of strong lines and the 
overall thermal shape, with peak near 1.65/xm, argue against the presence of an AGN in the 
T-Galaxy core. 
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Fig. 4. — Foreground corrected T-Galaxy polarization images over the central 22 x 22 arcsec 
region, (a. Top left) Contours of polarized intensity PI overlaid on a false color represen- 
tation of the Stokes / surface brightness excess [in mag arcsec"^] above sky, both stepped 
logarithmically. PI contour colors are referenced to the color legend at lower right in the 
panel; Stokes I false colors reference the upper left color legend, (b. Top right) Foreground 
and Ricean corrected percentage polarization Pcorj with contours linearly stepped from zero 
through 3%, and masked beyond dp of 4%. (c. Bottom left) Foreground corrected polar- 
ization position angle P.A.cor? with contours linearly stepped from +90° to +210°. (d. 
Bottom right) Stokes / (false color image) with white polarization vector overlays encoding 
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Fig. 5. — The T-Galaxy, as seen in 2MASS {H — K) colors after foreground correction, (a. 
Top Left) Color Excess E{H — K) map. Note the symmetric arm-like regions of redder 
excess colors offset to the East and West of center, (b. Top Right) Color Excess uncertainty, 
from 0.002 mag at center to 0.142 at the edge. (c. Bottom Left) Deprojected and rotated 
color excess map, showing reddened center, arms and bluer regions outside the arms. The 
axes are labeled for linear offset at the T-Galaxy distance, (d. Bottom Right) T-Galaxy 
zone classifications, based on colors, color uncertainties, and polarization uncertainties, as 
described in Section 5.1 and listed in Table [H Zone designations, from center outward, are 
the 'Core', 'Bulge', 'Interarm', 'Arms', and 'Outer'. 
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Fig. 6. — Comparison of percentage polarization in H band versus color excess in {H — K). 
Colored, crossed error bars locate the T-Galaxy foreground-corrected zone averages and la 
dispersions for the five zones defined in Figure |5ld. Black open diamonds are individual 
Taurus field stars with high polarization significance; black open triangles are color excess 
zonal averages of binned low polarization significance stars (see Section 5.3). Both sets of 
stars have been reddened and polarized by the Taurus Molecular Cloud comple x. The solid 



black curve labeled 'JKD - Milky Way' is the approximate behavior found by iJones et al. 



( I1992I ) toward Galactic stars while the gray band surrounding it captures la deviations of 



the stars defining the curve. The enclosed regions labeled 'NGC 4565' and ' NGC 891' en - 



compass points mapped in these galaxies into the Pr versus E{H — K) plane by 



JonesI f |l997h . 



(A color version of this figure is available in the online journal.) 
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Fig. 7. — Comparison of the zone-averaged percentage polarizations and color excesses, 
plotted versus zone in the T-Galaxy. Polarizations are shown as open diamonds, color 
excesses as open triangles. Polarizations generally increase from the Core to the Outer zone. 
Vertical error bars represent la uncertainties of the means. 
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Fig. 8. — Comparison of the polarization position angle, and numbers of pixels, versus the 
position angle for each of 25 azimuthal wedges in the T-Galaxy P.A. cor image, measured 
from the photometric center. Horizontal axis is position angle to the center of each azimuth 
wedge of pixels, measured E from N, offset by 270°. Left vertical axis is pixel averaged 
polarization P.A. Right vertical axis is number of pixels in each bin, shown as the lower 
curve of connected black triangles. T-Galaxy major axis is oriented along 142°, seen as the 
increases in pixel numbers near 140° and 320°. The polarization P.A. means are plotted 
for each zone (as colors) and for all the zones together (black diamonds). Core pixels are 
represented as a dashed violet line and error bar. Two gray diagonal bands show behavior 
for centrosymmetric scattering and for a radial pattern. 



(A color version of this figure is available in the online journal.) 
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Fig. 9. — Extinction and reddening of stars in the T-Galaxy region, as revealed using the 
RJCE method, (a. Left Panel) Color-color diagram for 358 stars located within 10 arcmin 
of the T-Galaxy an d having JHKL'M' m agnitudes. Horizontal {H — M') color is used in 



the RJCE method (IMajewski et al. 



201 ll ) for deter mining extinctions. Ve rtical axis shows 



he 



Bessell fc Bret 



Indebetouw et al 



( 19881) loci of colors for 



( I2OO5I ) reddening vector 



(J — K) color. The solid and dashed curves are t^ 
unextincted giants and dwarfs, respectively. The 
is shown as the angled, dot-dashed line, with A^J values (in mag) shown. The stars in the 
field surrounding the T-Galaxy exhibit significant extinction, (b. Right Panel) Histogram 
of RJCE extinctions determined from the stars shown in the left panel. Note the lack of 
any second strong extinction population that would indicate a second extincting layer. The 
direction to the T-Galaxy appears to pass through extinction associated only with the Taurus 
Molecular Cloud complex at about 150 pc distance. 
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Table 1. Polarization and Color Excess Properties by Region 



Region 


N 


CTp 


''■e(h-k) 


E{H - K) 


<P> 


< E{H -K)> 




(pixels) 


(%) 


(mag) 


(mag) 


(%) 


(mag) 


Core 


4 


< 0.12 


< 0.02 




0.83 (0.25) 


0.215 (0.017) 


Bulge 


26 


0.12 - 0.2 


< 0.03 




0.96 (0.21) 


0.192 (0.011) 


Interarm 


82 


0.2 - 0.7 


< 0.10 


< 0.17 


1.26 (0.14) 


0.110 (0.011) 


Arms 


61 






> 0.17 


1.10 (0.15) 


0.229 (0.011) 


Outer 


123 


0.7 - 1.2 


< 0.20 




1.08 (0.14) 


0.135 (0.015) 
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Table 2. Extinction Estimates toward the T-Galaxy 



Method 


Quantity 


Value 


Equiv. Av 
(mag) 


E{H - K) 
(Mag) 


Ref.'' 


IRAS/COBE 


Ah 


0.676 mag 


4.127 


0.260 


SFll 


IRAS/COBE 


Ah 


0.867 mag 


4.989 


0.314 


SFD98 


CO/"CO 


Nh, 


2 - 2.5x10^1 cm-2 


2 - 2.5 


0.13 - 0.16 


GHN08 


NICE 


Av 


1.80±0.17 mag 


1.80 


0.113±0.011 




NICER 


Ak 


0.10 - 0.15 mag 


1.8 - 2.8 


0.11 - 0.18 


LLAIO 


NICER/3 


Av 


2.18±0.08 mag 


2.18 


0.137±0.005 




NICER/5 


Av 


1.97±0.09 mag 


1.97 


0.124±0.006 




RJCE 


Av 


1.96±0.08 mag 


1.85 


0.123±0.005 




Adopted 






2.00±0.10 


0.125±0.009 





''SFll; Schlafly & Finkbeiner 2011 (from NED); SFD98: Schlegel, Finkbeiner & Davis 
1998 (from NED); GHN08: Goldsmith et al. (2008) - our visual interpolation; LCAIO: 
Lombajrdi, Lada, & Alves - our visual interpolation 



